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Comment on "Monomer Dynamics in 
and Single-Stranded DNA Polymers" 
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Despite the success of the bead-spring models p) in 
the description of the dynamics of polymers some prob- 
lems remain unsolved. Recently, a new "puzzle" has been 
reported 0. In the motion of individual monomers 
within polymer coils, single- and double-stranded (ss and 
ds) DNA, has been observed using fluorescence correla- 
tion spectroscopy. The mean square displacement (MSD) 
of the end monomer was interpreted within the Rouse (R) 
and Zimm (Z) models for flexible polymers. According to 
Q, an agreement with the R model has been found for 
dsDNA, while ssDNA followed the expected Z-type ki- 
netics. It is discussed in this Comment that the kinetics 
observed for long dsDNA is far from the R one and the 
macromolecule behaves rather as a Z polymer. The new 
puzzle is only due to an improper use of the RZ theory. 

The quantity indicating whether the dynamics is of the 
Z or R type is the draining parameter. For the polymer 
internal modes it can be defined as h(j>) = t p r/t p z oc 
p -1 / 2 , where p is the mode number, for the diffusion 
of the coil it is given by Dz/Dr. Here t v i are the re- 
laxation times in the corresponding model and Di the 
diffusion coefficients of the coil 0. The "pure" Z and 
R models correspond to infinitely large and small h(l), 
respectively. For the DNA parameters used in Q, h(l) 
is neither small nor large. So, for dsDNA with 23100 
base pairs (bp) h(l) w 2 and Dz/Dr k 4, so that the 
R dynamics cannot be expected. The experiments thus 
should be analyzed with no preliminary assumption on 
a specific dynamics. In Q the data were compared to 
the MSD given by the t 1 / 2 and t 2 / 3 laws for the R and Z 
polymers [l|,|3j. However, these laws assume continuous 
distribution of the internal modes and are only rough ap- 
proximations (with incorrect time dependence for both 
long and short times) to the more general expression, 
which for the end bead reads 


(r 2 (t)} 


6Dt + 4:N(a/ir) 2 ^2p- 


[1 - expH/T p ] , (1) 


"R modes" (with h(p) < 1, which for the experiment 
means p > 4), contribute only 13% to (r 2 (oo)) due to 
the internal modes. In the interval 0.02 - 10 ms these 
modes produce from 95 to 34% of the total internal part 
of MSD. Moreover, the diffusion contribution cannot 
be neglected as in Ref. [2]. For example, for 23100 bp 
DNA in the interval 10 - 200 ms (where the Z behavior 
of dsDNA has been identified |2|) the diffusion term 
contributes from 20 to 60% to the total MSD. 

Numerical calculations using Eq. (1) show that for 
the parameters from Q the correspondence between 
the data [2] and the theory is unsatisfactory: neither 
the R nor the Z model corresponds to the more general 
RZ model. The polymer parameters thus should be 
different. We have obtained them by fitting the RZ 
model (1) to the data 0- The diffusion term was kept in 
the consideration. The fit was performed with the con- 
ditions Na < L and b < a/2. The contour length L was 
calculated using the known distance between the base 
pairs along the chain (0.34 nm). Fixing N — 78, which 
corresponds to the same distance (100 nm) between 
the beads as in Q, we obtained a = 99.1 nm, b = 49.5 
nm and h(l) ~ 8.6 that indicates the (approximate) 
applicability of the Z model. When the R limit of Eq. (1) 
was fitted to the data, the obtained parameters yielded 
the Kuhn length (defined as I = Na 2 /L), which was 
only 30 nm while the accepted value is about 100 nm, 
in agreement with the RZ case (I « 98 nm). This is one 
more argument against the validity of the R model for 
the dsDNA dynamics. In summary, as distinct from the 
conclusion in 0], the RZ (being predominantly Zimm) 
model should be preferred. An analogous consideration 
of the ssDNA data [2] showed that also this polymer 
follows essentially the Zimm dynamics (with hCVj ~ 20 
for the 6700 bases ssDNA), in accordance with [2j. 
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ing us with the experimental data Q. This work was 
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D R + Dz is the Kirkwood diffusion coefficient PACS numbers: 82.35.Lr, 82.35.Pq, 82.37.-j, 87.15.He 


'R 

and the relaxation rates are r" 1 = r~„ + r^ 1 


■p ~ ' pR ~ ' pZ 
is the number of beads and a the mean square 

distance between the neighboring beads along the 

chain). Equation (1) is valid for the times t >> R 2 q/t] 

[4] , which is satisfied in [2] (p is the solvent density, 77 its 

viscosity, and R the hydrodynamic radius of the coil). 

Within the interval from 0.02 to 10 ms (where according 

to the R behavior was observed for dsDNA) the 

approximation (r 2 (t)) w 6D R t + (8k B Ta 2 t/ir 2 r/b) 1/2 (b 
is the bead radius) is inapplicable: for the parameters 
N = 68, a = 100 nm, b = 12.8 nm 2] the sum in Eq. (1) 
constitutes from 8 to 80% of its t x l 2 approximation. The 
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